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Abstract —This paper considers how to accurately estimate the 
minimum sum-rate so as to reduce the complexity of solving 
cooperative data exchange (CDE) problems. The CDE system 
contains a number of geographically close clients who send 
packets to help the others recover an entire packet set. The 
minimum sum-rate is the minimum value of total number of 
transmissions that achieves universal recovery (the situation when 
all the clients recover the whole packet set). Based on a necessary 
and sufficient condition for a supermodular base polyhedron to 
be nonempty, we show that the minimum sum-rate for a CDE 
system can be determined by a maximization over all possible 
partitions of the client set. Due to the high complexity of solving 
this maximization problem, we propose a deterministic algorithm 
to approximate a lower bound on the minimum sum-rate. We 
show by experiments that this lower bound is much tighter than 
those lower bounds derived in the existing literature. We also 
show that the deterministic algorithm prevents from repetitively 
running the existing algorithms for solving CDE problems so that 
the overall complexity can be reduced accordingly. 

I. INTRODUCTION 

Due to the growing amount of data exchange over wireless 
networks and increasing number of mobile clients, the base- 
station-to-peer (B2P) links are severely overloaded. It is called 
the ‘last mile’ bottleneck problem in wireless transmissions. 
Cooperative peer-to-peer (P2P) communications is proposed 
for solving this problem. The idea is to allow mobile clients 
to exchange information with each other through P2P links in¬ 
stead of solely relying on the B2P transmissions. If the clients 
are geographically close to each other, the P2P transmissions 
could be more reliable and faster than B2P ones. 

Consider the situation when a base station wants to deliver a 
set of packets to a group of clients. Denote V = {pi,..., p^} 
the packet set and K. = {1 ,,K} the client set. Due to the 
fading effects of wireless channels, after several broadcasts via 
B2P links, there may still exist some clients that do not obtain 
all the packets. Fig.Q]shows an example when the base station 
wants to disseminate 6 packets to 3 clients. In this figure, 'H t 
is the has-set that denotes the packets obtained by client j 
after B2P transmissions. Since the clients’ knowledge of the 
packet set may be complementary to each other, instead of 
relying on retransmissions from the base station, the clients 
can broadcast combinations of the packets they know via P2P 
wireless links so as to help the others recover the missing 
packets. We call Fig. |T| cooperative data exchange (CDE) 
system. For this kind of systems, there is a so-called CDE 
problem: how to find an efficient transmission strategy that 



{pi, P2. P3. P4, P5) {Pl,P2,P6} 

Fig. 1 . An example of CDE system: There are three clients that want to 
obtain six packets. The has-sets are "Hi = {pi, P2, P3, P4, Ps}, 2 = 
(pi 1 P2> P 6 } and H3 = (p 3 , P 4 > P 6 }■ 

achieves the universal recovery (the situation when all clients 
recover the entire packet set). 

Let r = (n,... ,rx) be a transmission strategy, where ry 
denotes the total number of linear combinations transmitted by 
client j. The CDE problem is usually expressed in the general 
form of |ffl| 

min/(r) 

r 

s.t. x u > n «ii.v*cjc, 

jex j&K\X 

E r i =a > (1) 

j'e/c 

where a denotes the transmission budget, the designated total 
number of transmissions among the clients. In <|T]), the expres¬ 
sion of /(r) is determined based on what kind of strategy is 
considered to be efficient. For example, if /(r) = w T r, the 
most efficient strategy is the one that minimizes the weighted- 
sum of transmissions HI; if /(r) = Yyjetc r i lo g(u)’ the most 
efficient strategy is the one that distribute the transmission rate 
in the fairest way (2). For most of the algorithms that solve 
a particular CDE problem, the value of a is assumed to be 
known a priori. For example, the algorithms proposed in fT|, 
El find minimum weighted sum-rate and fairest strategies, 
respectively, for a given value of a. 

However, problem (0} does not have solutions for all values 
of a. In fact, there exists a minimum sum-rate a* so that the 
constraint set in (|T|) is nonempty only if a > a*. But, the value 
of a* is usually unknown in advance. On the other hand, the 
algorithms proposed in (TJ, El can check the feasibility of 
a given value of a (i.e., whether there exists a strategy that 
can achieve universal recovery under budget a). Therefore, 
for the minimization problem over the minimum sum-rate 
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strategy set (problem (Q]) when a = a*), one can start with 
any arbitrary value of a and adjust it to the feasible values 
accordingly, e.g., setting a as the lower bound on a* derived 
in a, a, run the algorithms in a, a and increase the 
value of a until a feasible solution is found. By doing so, the 
algorithms in Q], @ should be called for many times until 
the final solution is found. But, the complexity of running 
algorithms proposed in fT[, J2) is not low. These algorithms 
involve submodular function minimizations, the complexity of 
which is at least 0(K 5, y + A" 6 )[] For this reason, estimating 
a value of minimum sum-rate that is close to a* prevents 
from repetitively running algorithms proposed in a, a 
and, therefore, helps reduce the overall complexity of solving 
problem ijT}- 

In this paper, we study how to accurately estimate the 
minimum sum-rate a* for both CDE systems that allow 
packet splitting (PS-CDE) and CDE systems that do not 
allow packet splitting (NPS-CDE). We start the study with an 
existing result derived in a, a: the crossing submodularity 
of | rijeACVr |. By using a necessary and sufficient condition 
given in ] 6 ] for the constraint set in O to be nonempty, 
we show that the exact value of a* can be determined by 
a maximization over all possible partitions of the client set /C, 
which is NP-hard Q. Instead of calculating the exact value 
of a*, we propose a deterministic algorithm to estimate a 
lower bound on a*. We show by experiments that the lower 
bound found by the deterministic algorithm is much tighter 
than the ones derived in 0, (4). In fact, the experiment results 
show that the deterministic algorithm returns the exact value 
of minimum sum-rate in most cases. In this paper, we also 
discuss how our results are related to those in |[T|, {5], |4)], [ 8 ) 
and prove that any PS-CDE problem can be converted to an 
NPS-CDE problem. 

II. System Model and Problem Statement 

Let V = {pi,...,pi} be the packet set containing L 
linearly independent packets. Each packet p, belongs to a 
field ¥ q . The system contains K geographically close clients. 
Define the client set as K, = {1,..., K}. Each client j G K 
initially obtains TLj C V. Here, TLj is called the has-set 
of client j. We also denote TLj = V \ TLj as the packet 
set that is missing at client j. The clients are assumed to 
collectively know the the packet set, i.e., Uj^TLj = V. The 
P2P wireless links between clients are error-free, i.e., any 
information broadcast by client j can be heard losslessly by 
client j' for all j' G K, \ {j}. The clients broadcast linear 
combinations of the packets in their has-sets in order to help 
each other recover the entire packet set V. For example, in 
the CDE system in Fig. |T] client 1 broadcasting pi + P 3 
helps client 2 recover P 3 and client 3 recover pi, and client 
2 broadcasting pi + p 6 helps client 1 recover p 6 and client 3 
recover pi. 

For a transmission strategy r = (ri,..., rx ), where rj 
denotes the total number of linear combinations transmitted by 

'There are many algorithms proposed for submodular function minimiza¬ 
tion problem. To our knowledge, the algorithm proposed in {2 has the 
lowest complexity 0 (K 5 "/ + K 6 ), where 7 is the complexity of evaluating 
a submodular function. 


client j, we call YhjeK. r o t ^ le sum-rate of strategy r. Denote 
lZ a the set that contains all transmission strategies that achieve 
universal recovery and have sum-rate equal to constant a. 7 Z a 
is nonempty only if a is greater or equal to the minimum 
sum-rate a*. For example, consider the CDE system in Fig. Q] 
and assume it is an NPS-CDE system. Since r G , a G No, 
i.e., we only consider nonnegative integer values of a. In this 
system, a < 3, lZ a = 0, i.e., there is no strategy with a 
sum-rate less than or equal to 3 that can achieve universal 
recovery. But, if a = 4, there exist three strategies that achieve 
universal recovery: IZ 4 = {(3, 0,1), (3,1, 0), (2,1,1)}. For 
example, in strategy (3,0,1), client 1 broadcasting pi + P 3 , 
P 2 + P 4 and P 5 and client 3 broadcasting p 6 can achieve 
the universal recovery. Therefore, the minimum sum-rate is 
a* = 4. Assume Fig. |T] is a PS-CDE system. Since r G K^, 
a G K.+, i.e., we consider the nonnegative real values of a. It 
can be shown that a* = 3.5 with the corresponding strategy 
set being 1Z 3 . 5 = {(2.5, 0.5, 0.5)} and lZ a = 0 if a < 3.5. 
Note, the way to implement strategy (2.5, 0.5, 0.5) is to break 
each packet into two chunks, e.g., the has-set of client 3 would 
be {p^, Pg 2 ), p!, 1 '', p^, pg 1 "*, Pg 2 ' }• Letting client 1 transmit 
5 chunks and both client 2 and client 3 transmit 1 chunk is 
sufficient to achieve the universal recovery. In this case, the 
packet normalized strategy is (2.5,0.5,0.5). 

It is shown in [11 that the CDE problems can be formulated 
by dU in general. Note, (|TJ is an integer programming problem 
for NPS-CDE systems. It can be seen from the above examples 
that for any CDE system there exist a minimum sum-rate a* 
such that 7 Z a is nonempty only if a > a*, i.e., problem © 
does not have a solution when a < a* since the constraint set 
is empty. Therefore, estimating the value of the minimum sum- 
rate a* helps determine whether an a is feasible or not, and, 
as discussed before, reduce the overall complexity of solving 
problem ©. 


III. Minimum Sum-rate 

As discussed in Section HI] estimating the minimum sum- 
rate a* is related to the nonemptyness of 7 Z a . In this section, 
we show a method to determine a* by studying the condition 
for the nonemptyness of 1Z a . We use a necessary and sufficient 
condition given in |j6) for lZ a to be nonempty to show how 
to determine a*. 


A. Constant Sum-rate Set 

Let r(X) = r .i anc * d e fi ne 


9(X) = 



( 2 ) 


3&K.\X 

It is proved in m that a transmission strategy r can achieve 
universal recovery if 


v(X) > g{X) (3) 

for all X such that X C 1C. The interpretation of © is: the 
information sent from X should be complement to that missing 
in K, \ X, i.e., the total number of packets transmitted by the 
clients in any set X should be greater or equal to the number 
of packets that are commonly missing at the remaining clients. 
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We can describe the constraint set in problem 0 as the 
constant sum-rate set 

lZ a = jr £ : v(X) > g{X),MX C fC, r (£) = a| (4) 

for a PS-CDE system. For an NPS-CDE system, r £ 
should be replaced with r £ N^ . Therefore, determining 
minimum sum-rate is equivalent to finding the smallest value 
of a such that lZ a is nonempty. 


B. Determining Minimum Sum-rate 

It is proved in Lemma 1 in (2) and Lemma 6 in Q] that 
g is a crossing supermodular functional Let ZC { 1 , 2 ,...} 
with 2 < \I\ < K and denote { X, },ax a partition of /C@ We 
can use the following proposition to determine whether TZ a is 
empty or not. 

Proposition 3.1 (nonemptyness oflZ a ): lZ a is nonempty if 
and only if 



(5) 

iex 


a -Yl (<* - 9{£ \ x i)) 

iex 

( 6 ) 


for all partitions {Xf\^x- 

Proof: This is a direct result of Theorem 2.6 in ©, which 
derives the necessary and sufficient condition for a crossing 
supermodular base polyhedron to be nenempty. The crossing 
supermodular base polyhedron in our case is TZ a . ■ 

Based on Pronosition l3.il we can determine the minimum 
sum-rate a* by the following theorem. 

Theorem 3.2: The minimum sum-rate is 

a* = max { g \x\ : a11 i X i}iGz\- ( 7 ) 

L iex ' ' 


for PS-CDE systems and 

g(>c\X) 
|Z|-1 


a = max 


iex 


: all partitions {Ti}j e xj 


for NPS-CDE systems, where \y \ is the ceiling function which 
returns the minimum integer value that is no less than y. 
Proof: Eq (0 in Proposition 13.11 can be rewritten as 


a 


y- g{K \ Xj) 

m-i 


(9) 


But, since g{X U^)> g{X) + g(y) for all 0 ^ X, y C /C, it 
can be shown that 

a«> 

iex ' ' iex 

So, lZ a is nonempty if and only if a > 'ffiex f° r 

all partitions. Therefore, for a PS-CDE system, the minimum 


2 A set function /: 2 K i-> R+ is supermodular if for all X,y C K. if 
f(X) + f(y) < f(X n T) + f(x u T) for all X,y C K such that 

x ny x -y fi,y - x Mid xe>y fK.®. 

3 A partition [X ] },xx of K, satisfies Xj f 0 , A',; f] X :j = 0 and U,cX'T: = 
/C for all i,j e X. 

4 Please see the proof in Appendix [Al 


sum-rate is determined by 0. For an NPS-CDE system, 
a* must be the smallest integer that is greater or equal to 
'Yh ei f° r partitions, which can be expressed by 

<nf ■ 

Consider what (0 and (0 mean in the CDE system. Assume 
that the clients can form groups, or coalitions, under the 
condition that any client can only appear in at most one 
coalition. Then, any form of coalition can be represented 
by {Xi}i e z, a partition of /C. Based on {Xi}i^x, consider 
determining the value of a such that lZ a is nonempty by using 
0. On one hand, r [Xf) > g{Xf) for all i £ I. Recall that a = 
r(/C) = X/iex r (Xi). We have the condition a > Siez 
which is exactly 0. On the other hand, r(/C\ Xf) > g(IC\Xi) 
for all i €2. We have Jfiei r (^ \ X *) - J2iex 9^ \ **)■ 11 
is equivalent to (|I| — l)a > which is exactly 

@. Since 0 should be satisfied for all subsets that are not 
equal to 1C, 0 and 0 should be satisfied for all partitions. 
But, according to (ITOt . 0 is more strict than 0. So, satisfying 
0 for all partitions is sufficient to determine a* as stated in 
Theorem 13.21 

Example 3.3: Consider the CDE system in Fig. |T| We have 
function g as 

5 (0) = 0,3({1}) = l,p({2}) = 0, ff ({3}) = 0, 

3({1,2}) = 3,3({1,3}) = 3,3({2 1 3}) = 1. 

By applying Theorem 13.21 we have a* = 3.5 for PS- 
CDE system and a* = 4 for NPS-CDE system. The cor¬ 
responding strategy sets are IZ 3.5 = {(2.5,0.5,0.5)} and 
IZ 4 = {(2,1,1), (3, 0,1), (3,1,0)}. It can be shown that lZ a 
is empty for all a < 3.5. 


C. Relationship with Existing Works 

In ®, it was shown that the minimum sum-rate is deter¬ 
mined by 


* r . / Yhiei ! Uj eXi | L 
a = L — mm <- 


m -1 

all partitions {<3x}i g zj@ 


( 11 ) 


One can show that it is exactly 0. Alternatively speaking, the 
authors in 0 derive the the minimum sum-rate for PS-CDE 
systems. 

For PS-CDE systems. Theorem 3 in 0 states that dividing 
packet into A '—1 chunks is sufficient to achieve the normalized 
minimum sum-rate with high probability. But, Theorem [30 in 
this paper establishes that dividing each packet into I\ — 1 
chunks is sufficient to achieve the normalized minimum sum- 
rate for sure. In addition, define 


g=(K-l)g(X). ( 12 ) 

We can also conclude from Theorem 13.21 that there exists a 
minimum integer valued a such that 

n a = {r £ r(X) > g(X),VX C /C,r(/C) = a} (13) 


5 It is found in ( 9 ) that the solutions of CDE and secrecy generation 
pro blem s can be found by solving the same linear programming (LP) problem. 
EqdJ was originally proposed in Go) to solve a secrecy generation problem. 
It is also a method to determine the minimum sum-rate of CDE systems (8). 
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Algorithm 1: Deterministic Algorithm 


Initiate the sum-rate lower bound: j3 *sr [ 

for k = 1 to K do 


Wr = {&}; 

for m = 2 to K — 1 do 


< = arg max{i?(/C\(VV m _iU{u}))—gr(W\{it}): 
IC\ Wm-i}; 

Wm = W m -1 U fuj,); 

/? = maxtf, 

endfor 
endfor 


- 1 ); 


Output /3; 


it G 


is nonempty, and that value is the non-normalized minimum 
sum-rate for the PS-CDE system. Therefore, any PS-CDE 
problem can be converted to NPS-CDE one (See the example 
below). For this reason, in the rest of this paper, we consider 
NPS-CDE problems only. 

Example 3.4: Consider the CDE system in Fig. [Q and 
assume it is a PS-CDE system. Let each packet be divided 
into K — 1 = 2 chunks. We replace g with g as 

g{%) = 0, p({l}) = 2,5({2}) = 0, p({3}) = 0, 
ff({l)2}) = 6 , <?({!,3}) = 6,p({2,3}) = 2 . 


By using © in Theorem 13.21 we get a* = 7. Therefore, 
the normalized minimum sum-rate is 7/2 = 3.5, which is in 
consistence with the result in Example 13.31 


IV. Tight Lower Bound on Minimum Sum-rate 

Problem © is equivalent to a minimum fc-partition problem. 
The NP-hardness of this problem is proved in 0, i.e., to 
directly determine the minimum sum-rate by Theorem 13.21 is 
intractable. Therefore, we propose a deterministic algorithm as 
shown in Algorithm 1 to approximate a*. In this algorithm, 
we initiate (3, the estimation of the minimum sum-rate, as 
the value of [ 1 bt t ^ le case -partition. We 

then update j3 by considering the value of [ J2iei 9 jip~i^ 1 
under a set of m-partitions where 2 < m < K — 1. From the 
following theorem, we show that the output 8 of Algorithm 1 
is a lower bound on a*. 

Theorem 4.1: For any CDE system, the output 8 of Algo¬ 
rithm 1 satisfies 8 < a*. 

Proof: We remark that W m in Algorithm 1 for all 
2 < to < K — 1 is generated based on Queyranne’s algorithm 
proposed in fm Here, W m and all j G K. \ W rn form a 
\K. \ Wm | + 1 -partition that contains \K \ W m | singletons. It 
is proved in fill that 


g(IC\W m )+g(JC\{jj) > g(fC\(W m \X))+g(lC\(XU{j})), 


for all j G K. \ W m , X C W m _i and 2 < m < K — 1. 
Therefore, 


g(IC\W m )+ J2 9(IC\{j}) 

j£lC\Wm 

>g(IC\(Wm\X))+g(JC\({f} UX)) 

+ E 9(lC\{j}), 


] is the maximum value of 


r s( ;c \ w ’")+^i«\w m 9(AL}) 

1 - e -> I |K\W m | 

f" g | i|E^ l over a su bset of \K \ W m | + 1 -partitions 

for all m. Therefore, based on Theorem 13.21 the output 8 of 
Algorithm 1 is a lower bound on minimum sum-rate a*. ■ 


A. Relationship with Existing Works 

In Lemma 2 in 0 , it is shown that the minimum sum-rate 
is lower bounded as 


a* > max |7f c | = maxg(1C \ {j}). (14) 

In Lemma 10 in 0 , it is shown that the minimum sum-rate 
is lower bounded as 


rE^I^-h 


rEj e K;s(£\{j})i 

I< - 1 


K - 1 


(15) 


But, it can be seen that the output 8 of Algorithm 1 satisfies 


8 > max { max { g (K 8 \ {j })+(/({/})}, 


[T,jeK9(£\ {i})l 
I< - 1 


}■ 


Since ma x j&K {g(lC \ {j}) + g({j})} > max ieK g(f C \ {j}), 
8 is tighter than the lower bounds given in 0, 0. To see 
how tight 8 is. we run the following experiments. 

Example 4.2: We vary the number of clients I\ from 3 
to 15 and the number of packets L from 6 to 30. For each 
combination of I\ and L, we repeat the procedure below for 
1000 times. 

• randomly generate the has-sets Hj for all j G K, subject 
to the condition U jejcHj = V', 

• get the minimum sum-rate a* by using the randomized 
algorithm proposed in 0; calculate the lower bounds 
given in 0 and 0 ; obtain the lower bound 8 by running 
Algorithm 1. 

We take the error as the absolute value of the difference 
between the lower bound and the minimum sum-rate a* and 
the average error as the mean error over 1000 repetitions. 
We show the average error of lower bounds in 0 , 0 and 
Algorithm 1 in Figs. [3 [3] and 0 respectively. It can be seen 
that the lower bound found by Algorithm 1 is much tighter 
than the ones in 0, 0 and, in most cases, Algorithm 1 finds 
the exact value of a*. 


B. Complexity 

It should be noted that running Algorithm 1 alone is not 
able to completely solve problem ©. To find a solution of 
problem © and also check if the universal recovery can be 
achieved with the sum-rate a = /3, one needs to run the 
algorithms proposed in 0 , 0 . But, running Algorithm 1 
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the number of clients K the number of packets L 


Fig. 2. The average error of lower bound in (3| when K varies from 3 to 
15 and L varies from 6 to 30. 



the number of clients K the number of packets L 


Fig. 3. The average error of lower bound in £0 when K varies from 3 to 
15 and L varies from 6 to 30. 


before those algorithms in ID, El helps reduce the overall 
complexity. 

As aforementioned, the complexity of algorithms in HI, 12 
largely depends on the complexity of solving an submodular 
function minimization problem. For example, the complexity 
of the algorithm in Q] is 0(K ■ SFM(K)), where SFM(K) 
is the complexity of running a submodular function mini¬ 
mization algorithm. Although it is proved that submodular 
function minimization problem can be solved in polynomial 
time, the lowest complexity of SFM{K) as proposed in 0 
is still 0 (AT 5 7 + K 6 ), where 7 is the complexity of evaluating 
function g. On the other hand, the algorithms proposed in Q), 
12 also require a certain value of a as an input. For example, 
the authors in fl] suggested use an initial value of a equal 
to the lower bound in |3), repeat running an algorithm with 
complexity 0(K ■ SFM(K)) and increase a until a feasible 
solution is found. Based on Fig. [3 it means that this algorithm 
should be called for approximately five times before it finds a 
feasible solution when K = 3 and L - 30. On the contrary, 
the complexity of of Algorithm 1 is just 0(K 3 7 )^ Based on 
Fig. [4] if we choose the initial value of a to be j3, the number 
of repetitions of running the algorithm in 12 , and hence the 
overall complexity, will be largely reduced. In fact, in most 
cases, the algorithm in ID will be called just once (since 
/3 = a* in most cases). In Fig. [3 we show the maximum 
error of /? collected in Example 14.21 It can be seen that the 
maximum error of j3 is one, which means that if the f3 ^ a*, 
the number of repetitions of running the algorithm in 0 ] will 
not exceed two. 


6 Algorithm 1 contains K(K — 2) iterations, and each iteration involves 
a maximization over client subset K \ W m —1 by evaluating function g. 
Therefore, the complexity is 0(K 3 7 ). 



the number of clients K the number of packets L 


Fig. 4. The average error of lower bound /I found by Algorithm 1 when K 
varies from 3 to 15 and L varies from 6 to 30. 



the number of clients K the number of packets L 


Fig. 5. The maximum error of lower bound /3 found by Algorithm 1 when 
K varies from 3 to 15 and L varies from 6 to 30. 


V. Conclusion 

In this paper, we showed that the minimum sum-rate a* in 
a CDE system could be determined by a maximization over 
all possible partitions of the client set. Instead of solving the 
maximization problem directly, we proposed a deterministic 
algorithm to estimate the lower bound on a* and showed by 
experiment that this lower bound was much tighter that those 
derived in the existing literature. The experiment results also 
showed that the output of Algorithm 1 was the exact value 
of a* in most cases and the maximum difference between 
it and a* is one, which could avoid repetitively running the 
existing algorithms for solving CDE problems so that the 
overall complexity could be reduced. 
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Appendices 


Appendix A 

Let such that X (~l y = 0. We have 

si(xuy)- g {X)- g {y) 


n 

- 

n «; 

- 

n 

je)C\(xuy) 


j€K\X 


jeic\y 


n n 

- 

n nf 

u 

n 



y) ieKVV 


n u ‘, 


n w;n(n«iun«i) 


je/c\(Wu:y) ie/c\(^uy) j£X jey 
> 0, 

i.e., g(Xljy) > g(,L)+g(V)[I|Then, for any partition {Xi}i G x 
of/C , 


g(/C \ Xi) > 9(Xi>)- ( 16 ) 




So. E ie z 9(JC \ *i) > m - 1) £ ieI </(*), i- e - 


E 

iez 


g(K\Xi) 

Pl-i 


>E^)- 

iGZ 


(17) 


7 Since g(0) = 0, g can be considered supermodular in all X such that 

m^x ClC. 

















